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1. Introduction 

Lattice studies of hadrons containing b quarks are important for several reasons. One major 
motivation is flavor physics, where non-perturbative calculations of hadronic matrix elements for 
electroweak transitions are required. Secondly, lattice QCD can predict masses of hadrons that 
have not yet been observed experimentally. A few singly-bottom baryons have been found so far, 
and more results are expected from the LHC. Most recently, the Q.}, baryon was discovered at 
Fermilab. There are now two incompatible results for its mass, obtained by the D0 ^ and CDF 
collaborations. Lattice QCD can contribute to resolve this discrepancy. 

A number of unquenched calculations of bottom baryon masses have been done recently [^, 
0, |6|, ^ (see also |j8|] for a review presented at this conference). It is important to perform 
independent determinations of the same quantities with different lattice formulations in order to 
test universality. In this work, the domain wall fermion action (with = 16, M5 = 1.8) is used for 
both the valence and sea u-, d- and s quarks, while the b quark is treated with non-relativistic QCD 
(NRQCD). Compai^ed to the static heavy-quark action, which was used in [|, H 0], NRQCD has 
the advantage that it is not limited to systems containing only a single b quark. Also, spin splittings 
which would vanish in the static limit can be calculated. 

This work makes use of the V = 24^ x 64 gauge configurations generated by the RBC and 
UKQCD collaborations [^. There are four different ensembles with pion masses ranging from 
about 672 to 331 MeV; the lattice spacing is approximately 0.1 Ifm. 



The form of the lattice NRQCD action used here is the same as in Ref. []10|], where the bot- 
tomonium spectrum was calculated on MILC gauge configurations with AsqTad sea quarks and 
Liischer-Weisz gluons. The RBC/UKQCD ensembles use different actions for both the sea quarks 
(domain wall) an the gluons (Iwasaki), and it is therefore a useful test of universality to compute 
the bottomonium spectrum again on these lattices before moving on to do heavy-light calculations. 



This was done in Ref. 1 1 1 ] . In addition to tests of the lattice actions, this work provided an accurate 
tuning of the bare b quark mass and an independent determination of the lattice spacing. The main 
results are summarized in Sec. ^below; the reader is referred to [ 1 1] for the details. 



Then, Sec. ^ goes on to describe the calculation of the bottom hadron spectrum, including B 
mesons, singly- and doubly-bottom baryons, and the triply-bottom Q.hbb- The heavy-light calcu- 
lations are still in progress, and here only results for amiight = 0.005, amstrange = 0.04 and limited 
statistics are shown. The full, chirally extrapolated results will be presented in a forthcoming pub- 
lication. 



2. Bottomonium 

The first step was the tuning of the bare b quark mass. When using NRQCD, all energies 
obtained from fits to hadronic two-point functions are shifted by some common constant, as the 
rest mass is not included in the theory. Thus, to tune the b quark mass it is convenient to consider 
the kinetic mass 

_ p^-[E{p)-E{0)f 

i[E{p)-Em ^^-'^ 

of the hadron. This is based on the relativistic continuum dispersion relation, which is in fact very 
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ami, 




T(25)-T(15) 
splitting 


ami 
0.005 


«25-i5 (GeV) 
1.740(25)(19) 


2.30 


4.988(12) 


0.3258(47) 


0.01 


1.722(38)(19) 


2.45 


5.281(13) 


0.3242(46) 


0.02 


1.708(92)(19) 


2.60 


5.575(13) 


0.3231(54) 


0.03 


1.72(12)(2) 



Table 1: Tji(15) kinetic mass and T(25) -T(15) 
splitting for three values of the bare b quark mass 
(lattice units). Errors are statistical/fitting only. 
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Figure 1: r\i,{lS) kinetic mass vs amy. Errors are 
statistical/fitting only. The Une shows a linear fit. 



Table 2: Results for the inverse lattice spacings of 
the different ensembles, obtained from the T(25) — 
T(15) splitting. 
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Figure 2: Radial and orbital energy spUttings in 
bottomonium. Errors are statistical/fitting only. 



close to the lattice dispersion relation in the case considered here: as demonstrated in Jl 1|], the 
speed of light is compatible with 1 within statistical errors of less than 0.25% for lattice momenta 
ap = n - In/L up to = 12; equivalently M\an shows no dependence on p within errors. 

Table |I| shows the kinetic mass of the rjh meson for three different values of amt,, on the 
ami = 0.005 ensemble. As can be seen in Fig. [l], the data are compatible with a linear depen- 
dence in the range considered. Fitting the function aMkin = A + B ■ ami, gives A = 0.489(25), 
B = 1.956(11). Also shown in Table is the T{2S) -T{\S) energy splitting, which is found to 
be nearly independent of ami,. The T(25') — T(15') splitting is furthermore expected to have very 
small systematic errors, and is therefore an ideal quantity to set the lattice scale by comparing to 



the experimental value of 0.56296(40) GeV []12[]. Then, using the experimental value of the rji. 



mass, 9.389(5) GeV [[13|], one can solve for the value of mi, that gives the correct kinetic mass in 
physical units. This gives 

ami, = 2.514(36). (2.2) 

Results for the lattice spacings of the four different ensembles, computed after tuning the b 
quark mass, are listed in Table ^. There, the first error given is statistical/fitting and the second is 
an estimate of the systematic errors (relativistic, radiative and discretization) due to the NRQCD 
action. 

Next, Fig. ^ shows results for the radial and orbital energy splittings. All masses have been 
determined by computing the energy difference to the T(15) and using the experimental T(15') 
mass as an input. The lattice scales were taken from Table |2[ Thus, the 2S and IS masses are not 
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Figure 3: Bottomonium S'-wave hyperfine splittings 
(energies relative to the T(15') and T(25') states, re- 
spectively). Errors are statistical/fitting only. 



Figure 4: Bottomonium P-wave spin splittings (en- 
ergies relative to the spin-average of the ;i;/ll^') 
states). Errors are statistical/fitting only. 



predictions here and hence no error bars are shown for them. The remaining energy splittings are 
in good agreement with the experimental results (lines). The sea quark mass dependence is found 
to be weak, as expected for sufficiently light quarks. 

Spin-dependent energy splittings were also computed and are shown in Figs. ^ and 0. Here, 
larger systematic errors are expected due to missing relativistic and radiative corrections as well as 
discretization errors. The latter are most severe for the (S-wave) hyperfine splitting, which is known 
to be sensitive to very short distances. The P-wave spin splittings shown in Fig. ^ are seen to be in 
relatively good agreement with experiment within the statistical errors. The \S hyperfine splitting 
was found to be 52.5 ± 1.5(stat) MeV on the most chiral ensemble, which has to be compared 
to the experimental value of 71.4+3 J (stat) ± 2.7(syst) MeV ||T3|]. In Ref. [|r|], the bottomonium 
spectrum was computed using a relativistic heavy-quark action on the same RBC/UKQCD gauge 
configurations. There, the hyperfine splitting was found to be only 23.7 it 3.7 (stat) MeV, a much 
larger deviation from experiment. 



3. Bottom mesons and baryons 

For the calculation of heavy-light meson and baryon masses, the set ofu/d and s valence quark 
domain wall propagators in use is an extension of the propagators that were computed and saved 
during the static-light calculation in [^. So far, only propagators on the ami = 0.005, anis = 0.04 
ensemble, with valence quark masses equal to the sea quark masses have been included in the 
ongoing NRQCD spectrum calculation. These quark masses correspond to pion and kaon masses 
of about 331 and 576 MeV, respectively. Note that also the strange quark mass is too large; the 
physical point corresponds to anis ~ 0.034 [^. 

The domain wall propagators have APE smeared sources. For the heavy-quark, NRQCD prop- 
agators are computed with both point and Gaussian smeared sources. Hadron correlation functions 
are then calculated for both point and smeared sinks and projected to zero momentum. They are 
fitted simultaneously in a fully correlated multi-exponential matrix fit, and errors are estimated us- 
ing bootstrap. The results shown below are from about 800 domain wall propagators. To increase 
statistics, correlators directed both forward and backward in time are computed. 
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Hadron 



Operator 



1+ 

2 

2 ' 2 

1 + 
2 

-7* 1+ 3" 

^hi 2 ' 2 

1+ 3- 
" '^'fo 2 ' 2 



D.h, D.I 



„ „^ J + 3 - 
O O* 1+3 



3 + 
2 



£abc iCYj)py 
£abc (C75)/3^ 
£abc {CYj)py 
£abc iCYj)py 
£abc iCYj)fiy 
£abc {CYj)py 
£abc iCYj)py 



q'; Q'a 
q'r' Q'a 

q'p Qa 
q'p Qa 

a h QC 

Q^pQ'^q'a 

QpQyQa 



Table 3: Operators for bottom baryons (C = 7472, 
nonrelativistic gamma matrix basis, m„ = m^). 



Splitting 



AM (MeV) AM (MeV) 
lattice experiment 



B*-B 

b:-b, 



^b ^b 
— D.h 



^bb ' 
O* - 

^hb 



■ ^bb 
■^bb 



48(9) 

49(4) 

25(25) 

18(16) 

19(10) 

24(14) 

38(9) 



45.78(35) 
46.1(1.5) 
21.2(2.0) 



Table 4: Heavy-light spin splittings in bottom 
mesons and baryons at ami = 0.005, am, = 0.04 
(preliminary; errors are statistical/fitting only) 



The structure of the baryon operators in use is shown in Table |3[ The b quark is denoted by 
Q, which is a 4-component spinor with vanishing lower components (in the nonrelativistic gamma 
matrix basis), since in NRQCD quarks and antiquarks are decoupled. The operators with Dirac 
matrix F = CYj have an overlap with both J = \ and J = \ states. At zero momentum, these 
contributions can be disentangled by multiplying the correlator with the projectors (5;y — ^YiYj) 
and 5777, respectively. 

As mentioned before, energies obtained from fits to correlators are shifted due to the use of 
NRQCD. Energy splittings are not affected. To compute the full hadron masses in a way that leads 
to only weak dependence on the bare b quark mass, the experimental value for the e.g. the T(15) 
or the B meson mass is used as an input parameter in the following way: 

M = £sim. + y (Mj,p. - (3-1) 
or M = £sim.+«fo(Mf^p.-£L.)> (3-2) 

where ri}, denotes the number of b quarks in the hadron, Ssim. is the simulation energy and M is the 



full hadron mass to be calculated. Results for the B meson masses, computed using (3. 1 ), are shown 
in Fig. ^ and the masses of singly- and doubly bottom baryons, for both methods (3.1) and (3.2), 
are shown in Figs. ^ and 0. Where available, the experimental values are indicated [12]; for the 
Q.h, both the D0 (black) and CDF (red) results are shown [[l], §]. Numerical results for various spin 
splittings are listed in Table 0; these are found to agree with experiment (where available) within 
the statistical errors. The hadron masses at the present values for the light quark masses tend to 
be slightly above the experimental results. Definitive conclusions can only be made after chiral 
extrapolation (and, eventually, after the inclusion of different lattice spacings and volumes). Note 
that the T mass shows little dependence on the sea quark masses, while the B has a light valence 
quark. Thus, ( p] ) and ( |3.2| ) lead to very different chiral behavior of M, which likely explains 
the discrepancies between the two methods seen at the present quark masses (for doubly bottom 
baryons. Fig. ^ the differences are enhanced since ri}, = 2). 
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Figure 5: B meson masses at ami = 0.005, ams 
0.04. Errors are statistical/fitting only. 



Figure 6: Singly bottom baryon masses at ami = 
0.005, ams = 0.04. Errors are statistical/fitting only. 



PRELIMINARY 



energy sliift from T ^- 
energy shift from B ^ 



>. 0.75 p 
9^ 0.5 h' 



0.25 



point - point 
point - smeared 

smeared - point • • 

smeared - smeared ^ — ' 



10 



time 



20 



30 



Figure 7: Doubly bottom baryon masses at ami = 
0.005, anis = 0.04. Errors are statistical/fitting only. 



Figure 8: £lhhb matrix correlator, effective energy 
plot (lattice units). 



The ^bbb baryon does not contain light valence quarks, and similarly to bottomonium, the 
dependence on the light sea quarks masses is expected to be weak once these are light enough. 
Thus, Eq. ( ^ ) is the better method for computing its absolute mass, and no chiral extrapolation is 
required. Also, since NRQCD is computationally cheap, one can go to very high statistics with little 
cost. An effective-energy plot for an Q-bbb matrix correlator from about 10^ NRQCD propagators 
on the ami = 0.005, aiUs = 0.04 ensemble is shown in Fig. ||. As can be seen, the signal is very 
good. The (unphysical) energy obtained from the fit is af'tiifci = 0.5527(12). Fitting an T correlator 
from the same propagators gives aEji^is) = 0.29786(20). Using the bootstrap method to properly 



take into account correlations, Eq. (3.1) then leads to 



Mn,,, = 14.3748(33) GeV (3.3) 

where the error is statistical only and includes the uncertainty in the lattice spacing (the latter was 
taken from Table |2|). The ^bbh mass has been estimated using various continuum methods, see 
[15, 16, 17, 18, |l^], and the production of the ^bbb athadron colliders has been studied in pO| , pT| ]. 
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4. Outlook 



The heavy-light calculations will be extended to include the other light quark masses, and 
chiral extrapolations will be performed. All calculations presented here are only for one lattice 
spacing, but the finer V = 32^ x 64 RBC/UKQCD gauge configurations will be included once they 
become available. This should allow more reUable estimates of discretization errors. 

With NRQCD, high statistical accuracy can be achieved for the O-bbh baryon, similarly to 
bottomonium. It should therefore be possible to study excited states also for the 0.hbb- 

Acknowledgments: Computations were performed at NCSA, NERSC, and Cambridge HPCS. 
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